“Calhoun 


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1993-06 


Adaptations to MICROPEP and ROCKET to 
allow performance evaluation of multiple 
grain and/or airbreathing motors 


McAtee, Aaron M. 


Monterey, California. Naval Postgraduate School 
http://hdl.handle.net/10945/39810 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 


: \§ D U DL EY research materials and institutional publications created by the NPS community. 
«iis eacica Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NNN KNOX appointed -- and published -- scholarly author. 


LIBRARY Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 


hittp:/fanwwenps.edu library Monterey, California USA 93943 





D-A268 798 RL 
CT ET ct (2) 


NAVAL POSTGRADUATE SCHOOL 
Monterey, California 





THESIS 


ADAPTATIONS TO "MICROPEP" AND "ROCKET" TO 
ALLOW PERFORMANCE EVALUATION OF MULTIPLE 
GRAIN AND/OR AIRBREATHING MOTORS 


by 


AARON M. MCATEE 


June 1993 





Thesis Advisor: D. W. NETZER 


Approved for public release; distribution is unlimited. 


nt san 





ONT aSSS 
Security Classification of this page 












REPORT DOCUMENTATION PAGE 
1b Restrictive Markings 


3 Distribution/Availability of Report 


2b Declassification/Downgrading Schedule Approved for public release; distribution is unlimited. 
4 Performing Organization Report Number(s) 5 Monitoring Organization Report Number(s) 

6a Name of Performing Organization 6b Office Symbol 7a Name of Monitoring Organization 

Naval Postgraduate School (if applicable) 31 Naval Postgraduate School 

6c Address (city, state, and ZIP code) To Address (city, state, and ZIP code) 

Monterey CA 93943-5000 Monterey CA 93943-5000 


6 Office Symbol 9 Procurement Instrument Identification Number 
(if applicable) 










































8a Name of Funding/Sponsoring Organization 
Naval Air Warfare Center, Weapons Division 










Address (city, state, and ZIP code) 
China Lake, CA 93555 


10 Source of Funding Numbers. 


Program Element No 





Work Unit Accession No 





11 Title (include security classification) ADAPTATIONS TO "MICROPEP” AND "ROCKET" TO ALLOW PERFORMANCE EVALUATION OF 
MULTIPLE GRAIN AND/OR AIRBREATHING MOTORS (UNCLASSIFIED) 


12 Personal Author(s) Aaron M. McAtee 
13a Type of Report 13b Time Covered 14 Date of Report (year, month, day) 15 Page Count 
Master’s Thesis From 8/92 To6/93} 1993,June 11 60 


18 Subject Terms 





























Di oa I 





Combustion, Performance, Codes, Computer, Rocket 





19 Abstract 
Adaptations to two existing rocket motor performance computer programs were made. MICROPEP, a FORTRAN program developed by the Naval Air 
Warfare Center Weapons Division, China Lake, California to evaluate theoretical performance of various propellants in rocket motors, was modified to allow 
calculation of the effects of non-ideal expansion and mixed shifting equilibrium-frozen composition nozzle flow on performance. In addition, the ability to 
handle vitiated air heaters, to calculate combustion chamber Mach number and to calculate stagnation pressure for airbreathing engines was incorporated. 
ROCKET, an internal ballistics program written by the Lockheed Corporation and modified by the Naval Air Warfare Center Weapons Division, was 
updated to purge the program of unused code and allow input of performance losses in both the combustor and nozzle flows. 

















| 
| 










20 Distribution/Availability of Abstract 21 Abstract Security Classification 

__ unctlassified/unlimited X same as report __ DTIC users Unclassified 

22a Name of Responsible Individual 22b Telephone (include Area Code) 22c¢ Office Symbol 
D. W. Netzer (408) 656-2980 AA/Nt 





Approved for public release; distribution is unlimited. 


Adaptations to "MICROPEP” and "ROCKET" to Allow Performance Evaluation of Multiple 
Grain and/or Airbreathing Motors 


by 
Aaron M. McAtee 
Lieutenant, United States Navy 


B.S.A.E., Auburn University 


Submitted in partial fulfillment 
of the requirements for the degree of 


MASTER OF SCIENCE IN AERONAUTICAL ENGINEERING 
from the 


NAVAL POSTGRADUATE SCHOOL 
June 1993 


Approved by: 





D. W. Netzer, Thesis Advisor 





DY J. Collins, Chairman 
Department of Aeronautics and Astronautics 





ABSTRACT 


Adaptations to two existing rocket motor performance 
computer programs were made. MICROPEP, a FORTRAN program 
developed by the Naval Air Warfare Center Weapons Division, 
China Lake, California to evaluate theoretical performance 
of various propellants in rocket motors, was modified to 
allow calculation of the effects of non-ideal expansion and 
mixed shifting equilibirium-frozen composition nozzle flow 
on performance. In addition, the ability to handle vitiated 
air heaters, to calculate combustion chamber Mach number and 
to calculate stagnation pressure for airbreathing engines 
was incorporated. 

ROCKET, an internal ballistics program written by the 
Lockheed Corporation and modified by the Naval Air Warfare 
Center Weapons Division, was updated to purge the program of 
unused code and allow input of performance losses in both 


the combustor and nozzle flows. 
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I. INTRODUCTION 


Many factors enter into the design of a rocket motor 
(required thrust, specific impulse, operation envelope, 
Susceptibility, safety, cost, flight time, weight, size). 
It is the job of the engineer to take the requirements and 
to choose a propellant and motor configuration that best 
meet these needs. In order to do this, predictions or 
calculation of performance must be made. Because of the 
high cost of test firing, it is desirable to use computer 
codes that will accurately predict the theoretical 
performance of various compounds, thereby eliminating the 
need to conduct test firings on all but a few acceptable 
candidates. In recent years these codes have become 
availebie in PC versions. 

After the propellants are selected, performance losses 
musi be estimated. With these losses, actual pressure-time 
and thrust-time profiles and specific impulse can be 
calculated for various solid propellant grain designs and or 
ramjet configurations. 

It was the intent of this investigation to modify two 
existing computer programs to cover a wider variety of 
performance and design problems. Specifically, MICROPEP 
[Ref. 1] (an equilibrium thermochemical "propellant 


evaluation code" developed by the Naval Air Warfare Center 





Weapons Division, China Lake, California) was to be modified 


to allow input of the desired nozzle area ratio and ambient 
pressure so that effects of non-ideal expansion could be 
accurately calculated. The performance value for shifting 
equilibrium flow to the nozzle throat and frozen composition 
flow afterwards is often a good approximation to the actual 
flow and was also to be incorporated. 

In order to incorporate airbreathing applications the 
inlet air temperature and composition, motor static 
pressure, nozzle contraction ratio and nozzle discharge 
coefficient are required inputs. With these inputs the 
program should calculate the vitiated air composition and 
heat of formation as well as the actual combustion chamber 
Mach number and stagnation pressure and the associated 
performance values. 

ROCKET is a solid propellant rocket motor internal 
ballistics code orginally developed by the Lockheed 
corporation. With igniter, propellant and insulation 
geometry and physical porperties input it calculates the 
pressure-time and thrust-time profiles. Some of the 
subroutines have never worked correctly and many unused 
portions of the code had never been removed. Documentation 
for the code is also very limited. Modifications to ROCKET 
were required to make the nozzle loss coefficient an input, 
to update the basic algorithm so that multiple grains could 


be burned at the same time and to incorporate plot routines 


for the pressure-time and thrust-time profiles. Average 


thrust, average pressure and total impulse were also desired 


output values. 


II. MICROPEP 


A. BACKGROUND 

The MICROPEP program is a PC based computer program 
developed by the Naval Air Warfare Center Weapons Division, 
China Lake, California and is used to estimate the 
performance of solid and liquid propellant rockets. It also 
has the capability of determining equilibrium and thermo- 
chemical properties at a specified pressure-temperature 
point. 

MICROPEP can be broken into three parts. The first part 
provides an evaluation of combustion product composition and 
chemical properties assuming constant pressure, equilibrium 
adiabatic combustion. The unmodified program allows input 
of up és ten ingredients and the output consists of the 
following combustion chamber properties: temperature, 
pressure, enthalpy, entropy, ratio of specific heats, 
chemical composition, molecular weight, heat content and 
speed of sound. 

The second part of the program provides an evaluation of 
the throat and exhaust compositions and chemical properties 
assuming isentropic, adiabatic expansion through the nozzle. 


The unmodified program calculates these properties for two 


separate conditions: 








1. shifting equilibrium flow - assumes that the 


reaction kinetics are fast enough to maintain chemical 
equilibrium throughout the expansion process 

2. frozen composition flow - assumes that the reaction 
rates are so slow that the exhaust composition is the same 
as the chamber composition. The output for the throat and 
exhaust are the same as for the chamber. 

The third part of the program provides an evaluation 
of nozzle performance parameters assuming one-dimensional 
flow, Mach 1 at the throat and exit pressure equal to 
ambient pressure. Output consists of the following for both 
frozen and shifting flow: specific impulse, ratio of 
specific heats for the chamber to throat process (the so- 
called "process gamma"), characteristic exhaust velocity, 
nozzle expansion ratio, density specific impulse, vacuum 
impulse, ratio of nozzle throat area to mass flow rate and 
nozzle thrust coefficient. An example of the output from 
the code is shown in Table 1. It should be noted that the 


units are a mix between English and SI. 


B. MODIFICATIONS 

None of the modifications made to the program used any 
new solution optimization algorithms. The existing 
subroutines were used, as needed, to affect the changes. A 
flow diagram of the MICROPEP program is shown in Figurel, 
with the modifications shown in bold. The changes made are 


as follows: 
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Table 1: 
NEWPEP SAMPLE OUTPUT 


#+e+ NEWPEP - Feb. 1990 ###¢# 


# xx * 05/16/93 * DH ** DENS ###* COMPOSITION *#éetes 

OXYGEN (GAS) 0 0.00001 20 

HYDROGEN (GASEOUS) 0 0.00001 2n 

INGREDIENT WRIGHTS (IN ORDER) AND TOTAL, WEIGHT (LAST TTlEM TN LIST) 
94.0000 6.0000 100.0000 


THE PROPRLLANT DENSITY IS 0.00001 LR/CU-IN OR 0.0003 GM/cCc 
NUMBER OF GRAM ATOMS OF EACH ELEMENT PRESENT IN INGREDIENTS 
5.952381 5.875000 0 


PR eHEE ENE EHEEERAR HEHEHE SEECHAMRER RESULTS FOLLOW tte thee eee eeednreeeaaene 


TCR) TCE) PAI) PC(ESI) ENTHALPY ENTROFY CE/CV SGAITIA RT/V 
3272. $430. 34.01 500.00 0.00 293.21 1.2N69 1.1361 7.284 


DAMPED AND UNDAMFED SPEED OF SOUND= 4061.859 AND 4061.860 FT/SEC 


SPECIFIC HEAT (MOLAR) OF GAS AND TOTAL} 11.590 11.590 
NUMBER MOLS GAS AND CONDENSFD- 4.669% 9.0000 


2.63700 H2O0 1.32936 02 0.45880 HO 0.11383 0 
8.93F-02 H2 3.77E-02 H 3.32F-03 HO2 6.49E-06 03 


THE MOLECULAR WEIGHT OF THE MIXTURE IS 21.417 


TOTAL HEAT CONTENT (298 REF) 71722.111 CAL/GH 
SENSIBLE HEAT CONTENT (298 REF) -1393.949 CAI,/GM 


POPPER H ERD EAS OSE EEEA EE EASAEEYHAUST RESULTS FOLLOWS tte teas bee e eat enaaaeeaes 
T(K) TCE) PATH) P(PSI) ENTHALPY ENTROPY CE/CV SGAMHA RT/V 
2192. 3486. 1.00 34.70 -f6.91 293.22 1.2923 1.1732 0.225 
DAMPED AND UNDAMPED SPRED OF SOUND- 3252.018 ANN 2252.018 FI/SFC 


SPECIFIC HEAT (MOLAR) OF GAS AND TOTAL- 11.349 ¥1.349 


NUMBER MOLS GAS AND CONDENSED- A.ARTOD n.nnno 
2.93A79 H20 1.43343 02 9.06255 HO 9.00648 9 
5.62F-03 #2 8.5NE-04 i 1.62F-04 HO2 


THE MOLECULAR WEIGHT OF THE MIXTURE 1S 22.482 
TOTAL HEAT CONTENT (299 REF) =1131.643 CAL/GN 
SENSTRLE NEAT CONTENT (298 RFEF)= 819.350 CAL/GN 





Table 1: (cont.) 
NEWPEP SAMPLE OUTPUT 


#444444444PERFORMANCE: FROZEN ON FIRST LINE, SHIFTING ON SECOND LINE*+##+#+# 


An exact method for determining throat conditions was used 

The frozen & shifting STATE gammas for the throat are: 1.2063 1.1360 
isentropic FXponent shown helow is the gamma for the chamber to throat 
PROCESS. 


IMPUhSE IS EX T* Pe cf Isp* OPT EX D-ISP AIM. EX T ADIT 
263.3 1.2089 2981. 19.83 5691.9 5.09 0.1 0.35390 1728. S84778. 
275.0 1.1360 3096. 19.63 5815.9 222.9 5.77 0.1 0.36161 2192. 770299. 


HEEHEEREREEAEH AEE EE EE EHHE EOE ATHROAT RESULTS FOLLOWS tte eee ade seed eseeetaetnrrse 


T(K) T(E) P(ATM) F(PSI) FNTHALPY ENTROFY CE/CV SGAMMA RT/V 
3096. 5114. 19.63 288.61 -16.14 293.21 1.2063 1.1360 4.250 


DAMPED AND UNDAMPED SFFED OF SOUND- 3929.313 AND 3929.314 FT/SEC 


SPECIFIC NEAT (MOLAR) OF GAS AND TOTAL= 11.621 11.621 
NUMBER MOLS GAS AND CONDENSEND- 4.6196 0.0000 

2.70089 H2O0 1.35026 O02 0.38104 HO 0.08796 O 
7.016-02 H2 2.71E-02 H 2.28E6-02 NO? 3.74F-06 03 


THE MOLECULAR WEIGHT OF THE MIXTURE IS 21.647 


TOTAL NEAT CONTENT (298 REF) -1618.754 CAL/GM 
SENSIBLE HEAT CONTENT (298 REF)-1297.611 CAL/GM 


#44* NEWPEP - Feb. 1990 ###3% 





NEWPEP 
Controlling Program 


Call Putin 


‘as vitiated 
air option 
chosen? 















Calf VITIATOR 
Input vitiated air Ingredients, 
desired pressure and 
temperature 




















Calculate composition and 
enthalpy of vitiated air at 
chamber inlet pressure and 
temperature 


















Save vitiated air data 
for input into 
combustor 


Input combustor 
ingredients and 
options 







Input chamber pressure, 
nozzle expansion ratio, 
ambient pressure and 

ingredient weights 






Was non-ideal 
expansion option 
chosen? 












Input chamber pressure, 
desired exit pressure 
and ingredient weights 

















Are finite combustor 
Mach number effects 
desired? 


Input nozzle 
contraction ratio 
(A4f/AS) 






Call EQUIVALENCE 
Catculate equivalence 
ratio of propellant 






Figure 1: 
PEP93 Flow Diagram 





Cali GUESS 
Come up with a crude guess 
of composition to get 


calculations off to a faster start 


Call HBAL 
Compute chamber composition, 
temperature, speed of sound, molecular 
weight, entropy, etc. for a given pressure 
and constant enthalpy 


Call OUT 
Print chamber conditions 


Call DESIGN 
Calculate chamber specific 
volume 
















Cail SBAL 
Calculate exhaust conditions 
‘for frozen composition flow 

and constant entropy 


Guess new exit 
pressure 


No 


Call DESIGN 


Calculate throat conditions 
assyming Mach 1 









Is calculated nozzle 
area ratio equal to 
desired area ratio? 







Y 
Compute performance values for on 


frazen composition flow 
Isp, C*, Cf, ete. 


° 








'as non-idea 
expansion option 
chosen? 


Figure 1: (cont.) 
PEP93 Flow Diagram 









Call SBAI_ 
Caiculate exhaust conditions 
for equilibrium flow and 
contant entropy 
















Call DESIGN 
Calculate throat conditions 
assuming Mach 1 


Compute performance 
values 


Was non-idea 
expansion 
gption chosen 





Guess new exit 
pressure 









Call 
MACHNUMBER 
Calculate chamber 
Mach number 














's calculated 
nozzle area ratio 
equal to desired 
area ratio? 





Yes 





No 






Ss this the fir 
time through the 
program? 





Calt OUT 
Print exhaust conditions 














Call SBAL 
Caliculate exhaust conditions for 
shifting equitibrium flow to the 
throat and frozen composition flow 
afterwards 


Call DESIGN 
Compute performance values 


as non-idea 
expansion option 
chosen? 













Yes 





Guess new exit 
pressure 















Was finite 
combustor Mach 
number option 
chosen? 









s cafculated 
nozzle area ratio — 
equal to desired 
area ratio? , 









Print all performance values 


Figure 1: (cont.) 
PEP93 Flow Diagram 
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1. The maximum allowable number of input ingredients 
has been increased to fifteen and the gas constant for the 
mixture is calculated. 

R = R/M 

2. The original program was extremely concise, leaving 
off units. Choosing clarity over brevity, units were added 
for all data. 

3. Performance calculations: 

a. Shifting equilibrium versus frozen composition 
flow- The theoretical performance of a rocket engine varies 
with the two assumptions. If the reaction rates of the 
chemicals are fast enough to maintain equilibrium throughout 
the nozzle then the additional energy gained by the 
reassociation of molecules and atoms will be reflected in an 
increased kinetic energy at the exhaust nozzle exit. 

Looking at the thrust and specific impulse equations 

Fo= mV ./G¢ - (De-Po) Ae 

Tap ba Fg, /mg,, 
it is apparent that theoretical performance would be highest 
for shifting flow and lowest for frozen flow. Since actual 
performance has been measured between these two conditions, 
modifications have been made to MICROPEP to calculate 
shifting flow to the throat and frozen flow afterwards. 
This accounts for the higher reaction rates due to higher 
temperatures and pressures coupled with a flow Mach number 


less than one prior to the throat. After passing through 


11 





the throat the flow accelerates past Mach 1 and the pressure 
drops rapidly to near ambient pressure, making it unlikely 
that the reaction kinetics will keep up with the flow 
velocity. Therefore, the assumption of frozen flow is 
reasonable. 

b. Non-Ideal Expansion - It can be easily shown 
that the maximum thrust occurs when the exit pressure is 
equal to ambient pressure [Ref. 2] 

F = mV./9, - (Pe-Po)Ae 
Taking the derivative of thrust with mass flow rate 
and ambient pressure constant gives 

GF = mdV,/g. + (Pe-P,) dA, + A.dP, 
Using the continuity equation and the momentum equation for 
ideal flow 

GF = (pAV) dV, + (Pe-Po) dA, + A,(-pVdV). 

GF = (Pe-Po) dAe 
therefore 

GF/dA, = 0 if De=P, 
In an actual rocket motor the exit area is constant, so the 
expansicn will only be ideal for one operating altitude. An 
option was added to input nozzle expansion ratio and the 
ambient pressure. Performance values for non-ideal expansion 


can now be calculated for different ambient pressures. 
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4. Exhaust kinetic energy 
One of the output values of NEWPEP was labeled 


adiabatic head. This value is calculated using the formula 


v./2 
and is therefore labeled kinetic energy of exhaust. 
5. Total and sensible heat content 
These were included in the original NEWPEP but the 


definitions were not clear. The present values are: 


a. Total Heat Content = Epn jhe; - Epnshes 
(without water condensation) 
b. Sensible Heat Content = Deni he, : Epnihes 
(with water condensed) 
6. Air-breathing motors 
a. Vitiated air 
The theoretical performance of the combustor 
varies with the composition and temperature of the incoming 
air. The vitiated air option requires the user to select 
the vitiator ingredients, inlet stagnation temperature and 
pressure. Using this information a pressure-temperature 
point calculation is made for the inlet composition and heat 
of formation. This "vitiated" air is then automatically 
entered as a user defined ingredient for the combustor 


calculations. 
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b. Combustor Mach number 

As combustor Mach number increases the 

difference between static and stagnation pressure increases 
De/P = (1+ ((y-1) /2) *M?) * (-y/y-1) 

The unmodified program assumes the Mach number is zero and, 
therefore, stagnation pressure is equal to static pressure. 

To account for finite Mach numbers, the nozzle 
contraction ratio can be input. If this option is chosen 
the program first calculates the nozzle and exhaust 
conditions for a chamber Mach number equal to zero. Then, 
using the calculated process gamma and given nozzle 
contraction ratio, the combustion chamber Mach number is 
solved for using an interval halving method. Using a lower 
bound of M = 0 and an upper bound of M = 1.001, the 
following equation is iterated until convergence is reached. 

O = M*(A,g/Ac) - [2/ (y+1) + (-y-1) / (y41) *M?) 1 * C (y+) /2 (7-1) 1 
The bisection method was chosen because of its guaranteed 
convergence and the small range of Mach numbers to be 
considered. 
c. Nozzle discharge coefficient 

One of the basic assumptions of MICROPEP is one- 
dimensional, isentropic flow. In order to more accurately 
predict mass flow rate through the nozzle, a streamline 
curvature correction (Cg) is now an allowable input. With 
the nozzle discharge coefficient the continuity equation 


becomes 
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Ms = DeALCgly/RT_] 9-3 12/ (y+) )* 0 (yt2) /2 (y-1)] 
= PrApcg/C* 
If an accurate value is not known the program default is 
1.0. 
7. Equivalence Ratio 
One definition of equivalence ratio is the ratio of 

oxidizer required for complete combustion of the fuel (to 
the highest oxidation state) to that available for 
combustion [Ref. 3]. Using this definition, subroutine 
EQUIVALENCE solves the following equation 

@ = [Eiey 2y x (Ng) o]/n, 
where K is number of reactant elements, 


nN: 


; is the number of moles of fuel element i, 


(ng), is the number of moles of oxygen required for 
complete combustion of one mole of fuel element i, 
no is the number of moles of oxygen available. 
Note: The program attempts to account for all common fuel 
species. See Appendix C: Subroutine EQUIVALENCE for a 


complete list of the elements incorporated. 


C. REVISED INPUT AND OUTPUT 
Appendix A contains the revised "users guide" for PEP93. 


A sample output is shown in Table 2. 
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Table 2: 
PEP93 SAMPLE OUTPUT 


#444 NEWPEP - Jan. 1993 *##*# 
4 xx * 
INGREDIENT MASS NE DENSITY COMPOS TT TON 
(gm) (kcal/kg) (kg/m**3) 
OXYGEN (GAS) 94.00 9 .0000 20 
HYDROGEN (GASEOUS) 6.00 -0 .0000 ral 


VOLUME PERCENT OF INGREDIFNTS (IN ORDER) 
94.0000 6.0000 


THE PROPELLANT DENSITY IS .00001 LB/CU-IN OR .0003 GM/CC 
THE EQUIVALENCE RATIO IS - 5066 


NUMRER OF GRAM ATOMS OF EACH ELEMENT PRESENT IN INGREDIENTS 
-§.952381 5.875000 O 


PPebRabaeed ERE REREE ERE EELEVACHAMBER RESULTS FOLSOWS Eee esse ede ass eases eaaeaans 


TEMP PRESSURE ENTHALPY ENTROPY CP/CV  SGAMMA ri/ns 
(R) (Mpa /ATtM/ PSI) (k3/kg) (kd /kg-®) (Mpa/kmol) 
3271.9 3.447/ 34.02/ 500.00 (00090 «12.268 «= 1.2069: 1.1362 73831.150 


DAMPED AND UNDAMPED SPRED OF SOUND- 1238.061 AND 1238.062 m/sec 


SPECIFIC HEAT (MOLAR) OF GAS AND TOTAL= 484.92 ABA.922 J/kg K 


NUMRER MOIS GAS AND CONDENSED= 4.6693 .onan 

(+-liquid, &=so0lid) 
2.63700 H20 1.32935 02 -45881 HO -11383 0 
@.93E-02 H2 3.77F-02 3.32F-03 HO2 6.A9E-06 03 


THE MOLECULAR WEJGHT OF THE MIXTURE IS 21.417 qm/mole 
THE GAS CONSTANT IS 388.22 J/kg-K 


TOTAL NEAT CONTENT (298 REF) - 7205.281 kJ/kg 
SENSIBLE HEAT CONTENT (298 REF) = §58632.259 kI/kg 


bee eee abe eee eee EHERHSSSSEERHAUST RESULTS FOLLOWS Ce eee sedeseesasdaasrseaas 


TEMP PRESSURE ENTHALPY ENTROPY CP/CV SGARMA Pi/ni 
(R) (MPa/ATM/FST) (kJ/kg) (kJ /kg-K) (MPa/kmol) 
2191.9 .101/ 1.00/ 14.70 -3636.2370 12.268 17.2123 1.1733 2278.684 


DAMPED AND UNDAMPED SPRED OF SOUND= 991.215 ANT 991.215 m/sec 


SPECIFIC NEAT (MOLAR) OF GAS AND TOTAL= ATA .82 ATA .A22 I/kqQ-K 
NUMBER MOLS GAS AND CONDENSED- 4.4479 .9000 
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Table 2: (cont.) 
PEP93 SAMPLE OUTPUT 


(*=liguid, &=solid) 
2.93879 H20 1.43343 02 .06255 HO .00647 0 
5.62F-03 H2 &.50E-04 H 1.62F-04 HO2 


THE MOLECULAR WEIGHT OF THE MIXTURE IS 22.483 gm/mole 
THE GAS CONSTANT IS 369.81 J/kg-K 


TOTAL HEAT CONTENT (298 REF) > 4734.795 kJ/kg 
SENSIBLE HEAT CONTENT (298 REF) = 3428.162 kJ/kg 


REFERER HEREEEEEERARERARERERHERHEAHEREEEREHEEHREEEEEEEREEEREEEEHEEEH EFA E ALAA EEAH EE 


An exact method for determining throat conditions was used 
The frozen & shifting STATE gammas for the throat are: 1.2063 1.1360 
GAMMA NU shown below is the gamma for the chamber to throat PROCESS. 


444404444 *¢PERFORMANCE: FROZEN ON FIRST LINE, SHIFTING ON SECOND LINE, ***+4¢ 
seese sees ees eSSHTETING TO THROAT/FROZEN AFTERWARDS ON THIRD LINR# Asse eeeeese 


SPECIFIC GAMMA T* re Ct TSP* Ae/A* D-ISP At/m Te Cf 
TMPULSE NU (gm-s/ (cm*##2/ 
(sec) (K) (MPa) (m/s) (sec) cm* #3) kg/s) (KR) 
263.3 1.2088 2962. 1.938 1733.3 5.099 2.017 $.0272 1728.1.4895 
275.0 1.1360 3096. 1.990 1772.7 222.9 5.768 2.107 5.1415 2192.1.5215 


267.2 1.1360 3096. 1.990 1772.7 222.9 5.110 2.047 5.1415 1816.1.4781 


SEREERHEEEEAEHHEREAEREEEESEPSEATHPOAT RESULTS FOLLOW Fete teee bese eee saanteaanat 


TEUP PRESSURE ENTHALPY ENTROPY CP/CV  SGAMMA PI/ni 
(RK) - (MPa/ATM/PSI) (kJ/kg) (kI/kg- Kk) (MEa/kmol) 
3096.5 1.990/ 19.64/ 288.61 -675.5911 12.268 1.2334 42.1441 2193.979 


DAMPED AND UNDAMPRD SPEFD OF SOUND= 1211.075 AND 1211.075 m/sec 


SPECIFIC HEAT (MOLAR) OF GAS AND TOTAL= 439.30 439.299 J/kg-K 
NUMBER MOLS GAS AND CONDENSFD- 4.6196 .oonn 
(*=ligquid, &-solid) 
2.70090 H20 1.35026 02 - 38103 HO -08795 0 
7.016-02 2 2.71F-02 2.2RF-0% HO? 2.74F-96 M21 


THE MOLECULAR WEIGHT OF THE MIXTURE IS 21.647 gm/mole 
THE GAS CONSTANT IS 384.09 J/kq-R 


TOTAL HEAT CONTENT (298 REF) = 6772.908 kJ/kg 
SENSIALE HEAT CONTENT (298 RFF) = 5429.259 kJ/kg 


FROZEN & SHIFTING KINETIC FNERGY OF EXHAUST 584773. 770295. m*#*2/s#42 
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III. ROCKET 


A. BACKGROUND 
1. ROCKET is a PC based program originally written by 
the Lockheed Corporation and then modified by the Naval Air 
Warfare Center Weapons Division. The intent of the original 
program was to calculate the pressure-time and thrust-time 
profiles of a multiple grain rocket motor using the 
following assumptions: 
a. Combustion products are ideal gases. 
b. Burning rate follows Veille’s Law 
r = ap,” 
c. Effects of mass addition and erosive burning 
are negligible. 
a. Inertia of the chamber gases is negligible. 
e. Variation of C* with pressure is described by 
Ce = C*pre¢ [Preg/1000]* 
f. Propellant burning rate temperature 
sensitivity is described by 
r= ryorexp (0.010, (T.-70) ] , where o, = (1-n)I 
g. Properties of the chamber gases are mass 
averages of the individual gases. 
2. Allowable inputs were: 


a. Number and type of grains 


(1) Igniter 





b. 


(2) 
(3) 


(4) 


Insulation 
Propellant 


Igniter Motor (for large rockets) 


Motor parameters 


(1) 
(2) 
(3) 
(4) 
(5) 
(6) 


(7) 
(8) 
(9) 


Nozzle throat and exit areas 

Total motor volume 

Ambient temperature and pressure 
Throat closure blowout pressure 
Initial motor pressure 

Nozzle discharge coefficient and half 
angle correction 

Ratio of specific heats 

Throat area design pressure 


Throat erosion rate(radial) 


Grain Parameters 


(1) 
(2) 
(3) 
(4) 
(5) 
(6) 
(7) 
(8) 
(9) 


(10) 


Burning rate 

Burning rate exponent 

Temperature sensitivity 

Reference pressure for C* correction 
Characteristic exhaust velocity 
Pressure correction exponent for C* 
Density 

Ignition time 

Flame spread ignition delay time 


Burning surface area 
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3. Using these inputs the program calculates chamber 
pressure by integrating the following mass conservation 
equation: 

dp./dt = (1/v9) [RT(ZppA,ap."-DAgnIc/C*) -Pe (dug/dt) } 
where v, is the instantaneous gas volume. Then the thrust 
coefficient is calculated using 
Cp = {2y7/ (y-1) * (2/ (yt1)) * LU (yt1) / (y-1)] 
*(1- (De/Pe) *(y-1)/y] }9°8X + (De-Pam) Ae/ (PeCaAen) 
With the values of pressure and thrust coefficient known, 


thrust was solved for using 


P = CppArnAcg 


B. MODIFICATIONS 

A flow diagram of the program is shown in Figure 2 anda 
modified user’s guide is given in Appendix B. 

1. Calculation of actual thrust and specific impulse. 
Several factors contribute to the reduction in theoretical 
thrust and specific impulse from the theoretical values 
(F’, Isp’). In general these can be put into two separate 
categories; losses due to incomplete combustion in the 
chamber and nozzle losses (two-phase flow, divergence, 
etc.). The ROCKET program has been modified to allow these 
losses to be accounted for using two efficiencies, yo. and 
Ic£- 

a. Combustion Loss Coefficient (n¢.) 
Nc» is primarily a function of the gas residence 
time and is determined empirically [Ref. 4]. Typical values 
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MAIN 
Controlling program 





Call INPUT 
Read in general motor 
description, grain parameters 


Caiculate propellant 
volume and weight 









Correct C* for 
actual chamber 
pressure 


Calculate 
throat area 





Ss input throa 


atea positive?, based on 


input 
pressure 


Yes 


Input throat area 
assumed correct and 
program only runs once 


Call EVES 
Integrate pressure equation over 
small time increments 


Call DIFEQ 
Caiculate throat erosion, change in 
pressure, new chamber pressure 
and mass flow rate 


No Call PRINT 
Print current pressure, change in 
pressure, free volume, propeliant weight 
consumed, insulation weight consumed 
- and throat area 











's Chambe 


pressure - Yes 
ambient pressure 
< 10 psia? 
Figure 2: 
ROCKET Flow Diagram 


21 





Call GETAVG 


Calculate average pressure 
over total burn time 


Print average pressure 
and maximum pressure 


Calt THRUST 
Calculate thrust, 
coefficient of thrust, 
specific impulse and total 
impulse 





Call GETAVG 
Calculate average thrust 
over total burn time 


Print average thrust 
and total impulse 


Calt PLOT 
Plot thrust vs. time and 








pressure vs. time 


for aluminized solid propellants are between .96 and .995 
for residence times greater than 10 seconds. For these 
propellants yo. drops off rapidly for residence times below 
10 seconds. 
b. Nozzle Loss Coefficient (1¢¢) 
This can be expressed by [Ref. 5]. 

Nef = 1- (€ givt €xint €b1+ Esubt ferost tp) /100 
where 

€aiy = half angle correction for conical nozzles 

Exin = COYrrection for shifting equilibrium vs. 


frozen composition flow 


€bh, = boundary layer loss 

Esub = COrrection for flow around submerged 
nozzles 

€eros= Losses due to nozzle throat erosion 

Erp = velocity and thermal lag of condensed 


species in chamber 

For use in ROCKET the input value of C* should be C* 
theoretical times no, . Allowable inputs of nozzle 
discharge coefficient and half-angle correction have been 
replaced by the nozzle loss coefficient, and thrust is 
calculated using 

F = MceCePcAth 
2. Plotting routines 
A plotting subroutine has been added which 


automatically takes the time and thrust values, determines 
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the maximum value of thrust and plots the thrust-time 
profile scaled to maximum thrust. The same is done for the 
pressure-time profile. In addition, time, pressure and 
thrust are printed to a separate data file. 

3. Total impulse is calculated using the following 


formula: 


where 

F is the average thrust 

t, is the total burn time. 

4. Improved efficiency 
Since the original program was written by Lockheed 

for use on a mainframe computer, ROCKET has undergone some 
modifications. These changes resulted in a program that was 
inefficient because several input variables were undefined 
in the weve manuals. In order to increase efficiency, more 
than 25 variables and over 300 lines of code associated with 
these variables were removed. The result is a program that 
still gives the output discussed above and does indeed work 


well for more than one grain. 


C. REVISED INPUT AND OUTPUT 
A sample input file is given in Appendix B with a rocket 
motor description shown in Figure 3. A sample output file 


is shown in Table 3. 
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Figure 3: 
ROCKET Description 
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Table 
ROCKET SAMPLE OUTPUT 


1Spercent 


Motor Ferformance Program IBM-FC Version 1.0 


INPUT DATA... 
Normal Input is finished 


General Configuration PARAMETERS are ... 





3: 


Throat AREA 8.6577 
Exit AREA 105.4500 
Expansion RATIO 12.1799 
Total Motor VOLUME 8042.5 
AMBIENT Temperature 70.0 
AMBIENT Pressure 10.1 
Closure BLOWOUT 35.0 
Pzero 14.7 
CF Efficiency 9606 
Gamma 1.1443 
Throat DESIGN Pressur 1030.0 
Throat FROSION .00n5 
MOTOR HAS 4 GRAINS 
DESCRIPTION OF GRAIN 1 
BURN RATE 5800 
BURN RATF EXP. 3500 
PI SUB K - 2000 
BURN REF. PRESS. 1030.0 
C STAR 5106.7 
C STAR EXP. -0450 
DENSITY 0623 
IGNITION TIME 0000 
DELTA IGN. TIME 0000 
WEB BURN AREA 
.000 744.150 
-381 759.750 
- 762 771.250 
1.143 778.650 
1.524 781.950 
1.905 781.150 
2.286 776.250 
2.666 767.250 
3.047 754.100 
3.428 736.900 
3.809 715.600 
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Table 3: (cont.) 


ROCKET SAMPLE OUTPUT 


DESCRIPTION OF GRAIN 


BURN RATE 


BURN RATE EXP. 


PI SUB K 


BURN REF. PRESS. 


C STAR 
C STAR FXP. 
DENSITY 


IGNITION TIME 
DELTA IGN. TIME 


WEB 
- 000 
391 
762 
143 
524 
905 


3.809 


2 


BURN AREA 


JAA. 
759. 
771 
778 
781 
781. 
776. 
767 
754 
736. 
715. 


DESCRIFTION OF GRAIN 


BURN RATE 


BURN RATE FXP. 


FI SUB K 


BURN REF. PRESS. 


C STAR 
C STAR FXP. 
DENSITY 


IGNITION TIME 


DELTA IGN. 


THIS GRAIN IS INSULATION WITH AN ARLATION RATE OF 


WER 


lad Ud GD AD NS tot 
“ 
a 
n 


TIME 


150 
750 


.250 
-650 
-950 


150 
250 


.250 
-100 


900 
600 
3 


BURN AREA 


1488. 
1519. 
1542. 
1557. 
1563. 
1562. 
1552. 
1534. 
1508. 
1473. 
1431. 


300 


5800 
3500 
- 2000 
1030.0 
5106.7 
-0450 
0623 
. 0000 
.9000 


- 5800 
3500 
- 2000 
1030.0 
5106.7 
.0450 
- 0623 
0000 
- 0000 
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95000 


Table 3: (cont.) 
ROCRET SAMPLE OUTPUT 


DESCRIPTION OF GRAIN 4 


BURN RATE 1.0000 
BURN RATE EXP. - 0000 
FI SUB K -0000 
BURN REF. PRESS. 1000.0 
C STAR 5106.7 
C STAR EXP. -0450 
DENSITY -0623 
IGNITION TIME -0000 
DELTA IGN. TIME -0000 
THIS GRAIN IS AN IGNITER 
TIME WDOT 
-0000 1.4000 
: -2500 1.9000 
TOTAL PROPELLANT WRIGHT= 362.478 IGNITER WEIGHT- AL 
1 
TIME EC Boor FREE VOL P WEIGHT WDOTP INS WT wnoT AT 
.000 14.70 22173.60 2217.59 -00 13.55 00 -00 8.65770 
-001 27.75 27381.44 2217.75 -O1 16.58 -00 .00 8.65770 
-001 38.26 27371.20 2217.88 -O1 18.39 -00 2.42 8.65770 
001 51.96 29634.08 2218.07 .02 20.31 -01 3.24 8.65770 
-002 66.64 31490.92 2218.27 -03 22.03 -0] 4.11 8.6577) 
-002 76.89 32543.34 2218.41 04 23.09 -01 4.71 8.65773 
.002 87.46 33467.93 2218.55 -05 24.09 -O1 5.32 8.65771 
003 98.30 34278.57 2218.70 06 25.04 -01 5.95 8.65773 
-003 109.38 34986.84 2218.86 -06 25.94 -01 6.59 8.65771 
-003 120.68 35602.51 2219.02 -07 26.80 .02 7.24 8.65772 
-004 132.16 36134.02 2219.18 08 27.62 02 7.90 8.65772 
-004 143.80 36588.74 2219.35 09 28.41 -02 8.56 8.65772 
004 155.57 36973.17 2219.52 -10 29.16 -02 9.23 8.65772 
-005 167.45 37293.06 2219.70 -il 29.89 02 9.90 5.65772 
-005 191.48 37759.36 2220.06 -13 31.26 02 71.25 8.65773 
- 006 215.74 38023.35 2220.43 15 32.53 .03 12.61 8.65773 
.007 240.11 38114.57 2220.82 -17 33.72 .03 13.97 8.65773 
-007 264.49 38057.92 2221.22 .19 34.83 -03 15.32 8.65774 
-008 288.80 37874.59 2221.63 22 35.88 .04 16.66 8.65774 
-009 312.95 37582.82 2222.05 -24 36.86 .04 17.99 8.65774 
-009 336.89 37198.40 2222.48 26 37.79 -04 19.30 8.65774 
-010 360.55 36735.09 2222.92 29 38.67 05 20.59 8.65775 
.010 383.89 36204.92 2223.37 31 39.50 05 21.86 8.65775 
-O11 406.88 35618.45 2223.83 34 40.28 -05 23.11 8.65775 
O12 451.65 34312.61 2224.77 .39 41.73 .06 25.53 8.65776 
O14 494.66 32879.32 2225.74 44 43.04 .06 27.85 8.65777 
015 535.78 31366.61 2226.73 -50 44.22 07 30.06 8.65777 
016 574.94 29811.89 2227.74 -56 45.30 -08 32.15 8.65778 
-018 612.09 28244.19 2228.78 -62 46.27 .08 34.13 8.65779 
-019 617.24 26685.84 2229.83 68 47.16 09 36.00 8.65779 
020 680.42 25153.92 2230.91 74 47.97 09 37.76 8.65780 
N21 711.65 29661.14 2232.00 . 80 AB.71 -10 39.42 8.6578) 
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Table 3: (cont.) 
ROCKET SAMPLE OUTPUT 


TIME EC PDOT FREE VOL PF WRIGHT WDOTP INS WT WDOT Ar 
023 FAL.OL1 22216.91 2233.10 - 86 49.39 -11 40.97 8.65781 
024 768.55 20827.85 2234.22 -92 50.01 -11 42.42 8.65782 
-025 794.35 19498.47 2235.34 .99 50.58 -12 43.78 8.65783 
026 818.49 18231.46 2236.48 1.05 §1.10 -12 4S.0S5 8.65783 
029 862.10 15889.49 2238.79 1.19 $2.02 -13 A7.34 8.65785 
-032 900.05 13800.72 2241.14 1.32 52.80 .15 49.33 8.65786 
-034 932.96 11954.15 2243.51 1.46 53.46 16 51.05 8.65787 
-037 961.43 10332.46 2245.91 1.59 54.02 .17 52.54 8.65789 
-039 986.03 8916.11 2248.33 1.73 54.49 18 53.82 8.65790 
-042 1007.24 7684.14 2250.77 1.87 54.90 .19 54.93 8.65791 
-044 1025.51 6616.43 2253.22 2.01 55.25 -21 55.88 8.65793 
-047 1041.23 $693.39 2255.69 2.16 55.54 .22 56.69 8.65794 
-050 1054.76 4897.45 2258.17 2.30 55.80 -23 57.40 8.65795 
-052 1066.40 4212.11 2260.66 2.44 56.02 .24 58.00 8.65797 
-055 1076.41 3623.16 2263.16 2.59 56.20 .25 58.52 8.65798 
-060 1092.43 2683.66 2268.17 2.87 56.51 28 59.35 8.65801 
-065 1104.32 1993.59 2273.20 3.16 56.74 -30 59.97 8.65803 
-073 1116.71 1288.38 2280.79 3.60 56.98 34 60.61 8.65807 
-083 1126.81 740.77 2290.94 4.19 57.20 -38 61.14 8.65813 
.106 1137.17 266.38 2313.86 5.51 57.46 -A49 61.68 8.65825 
-183° 1148.30 113.77 2390.75 9.94 57.95 -85 62.25 8.65865 
-249 1155.79 -2965.13 2457.94 13.81 56.44 1.16 62.64 8.65900 
-255 1140.96 -2095.53 2463.65 14.14 56.19 1.19 61.88 8.65903 
-263 1128.01 -1334.14 2471.24 14.57 55.98 1.22 61.21 8.65907 
-273° 1117.77 -724.26 2481.10 15.15 55.82 1.27 60.68 8.65912 
-J11 1107.56 -36.76 2518.33 17.29 55.71 1.45 60.15 8.65932 
-322 1107.41 2.71 2528.26 17.86 55.73 1.50 60.14 8.65937 
-532 1118.03 55.27 2732.80 29.61 56.30 2.49 60.70 8.66047 
-736 1128.33 42.24 2934.31 41.20 56.81 3.45 61.24 8.66154 
-982 1138.41 41.30 3178.18 55.22 57.32 4.62 61.77 8.66282 

1.228 1148.53 41.29 3424.17 69.37 57.83 5.80 62.31 8.66410 

1.576 1158.84 26.43 3775.79 89.60 58.34 7.48 62.85 8.66592 

2.006 1169.16 12.866 4213.94 114.82 58.86 9.56 63.40 8.66816 

2.539 1175.18 11.30 4759.95 146.24 59.18 12.15 63.74 8.67094 

2.600 1175.38 -1.67 4823.12 149.88 59.18 12.45 63.75 8.67126 

3.112 1173.38 -4.36 5349.12 180.16 59.09 14.94 63.66 8.67393 

3.645 1164.48 -19.20 5893.98 211.51 58.65 17.52 63.22 8.676712 

4.054 1153.28 -33.32 6309.70 235.44 $8.08 19.50 62.66 8.67885 

4.382 1142.12 -33.88 6639.00 254.38 57.54 21.08 62.09 8.68056 

4.648 1130.45 -48.58 6904.06 269.62 56.94 22.35 61.50 6.68195 

4.894 1118.50 -A4B.42 7146.19 283.54 56.35 23.51 60.88 8.68324 

5.099 1108.47 -55.57 7346.10 295.03 55.84 24.47 60.37 8.68430 

5.283 1097.31 -61.91 7524.33 305.28 55.28 25.33 59.80 8.68527 

5.447 1087.16 -61.81 7681.28 314.29 54.78 26.10 59.27 8.68612 

5.611 1077.07 -61.42 7836.83 323.23 54.28 26.85 58.75 8.68698 

5.774 1066.62 -70.73 7990.99 332.08 53.75 27.61 58.21 8.68784 

5.938 1054.58 -74.37 8143.54 340.83 $3.14 28.35 57.59 8.68869 

6.102 1042.41 -74.03 8294.40 349.49 52.54 29.09 56.96 8.68955 
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cc] 
~ 
a 
by | 


- 266 
-360 
-360 
-361 
361 
-362 
-363 
.365 
-366 
-368 
371 
.376 
-381 
- 386 
-391 
.397 
-402 
407 
412 
417 
422 
427 
432 
-438 
443 
.448 
453 
- 458 
-463 
473 
484 
494 
-509 


ANDAAAMAAMNRARAAAHAANAHHADHAAANAANHRAKTRAARAARARHN 


AVERAGE FRESSURE 
MAXINUM PRESSURE 


TIME 
. 000 
-001 
-001 
-O01 
-002 
-902 
002 
003 
.003 
.003 
- 004 
004 


eC 
1030. 
1019. 
1006. 
990. 
974 
958. 
928. 
698. 
869 
B15. 
763 
670. 
588. 
516. 
453 
398 
349. 


22. 


32 
22 
86 
63 


.65 


93 
25 
54 


-79 


03 


-71 


56 
77 
95 


-90 
-53 


92 


-24 
-76 
- 86 


FDOT FREE VOL 
-73.41 8443.60 
-25894.06 8527.95 
-25580.20 8527.95 
-25167.67 8527.95 
-24761.79 8527.95 
-24362.40 8527.95 
-23582.91 8527.95 
-22828.35 8527.95 
-22097.94 8527.95 
-20706.56 8527.95 
-19402.71 8527.95 
-17036.25 8527.95 
-14958.36 8527.95 
-13133.84 8527.95 
-11531.85 8527.95 
-10125.25 8527.95 
-8890.23 8527.95 
-7805.84 8527.95 
-6853.72 8527.95 
-6017.73 8527.95 
-5283.71 8527.95 
-4639.22 8527.95 
-4073.35 8527.95 
-3576.49 8527.95 
-3140.24 8527.95 
-2757.20 8527.95 
-2420.88 8527.95 
-2125.58 8527.95 
-1866.31 8527.95 
-1438.77 8527.95 
-1109.18 8527.95 
-855.08 8527.95 
-578.79 8527.95 
FROM -000 SEC. 
WAS 1175.41 FST 
FE F 
.00 oo 
.00 00 
5.18 273.95 
4.87 417.95 
4.63 577.29 
4.50 690.89 
4.39 809.64 
4.29 933.04 
4.20 1060.59 
4.11 1191.83 
4.04 1326.36 
3.97 1463.80 


Table 3: 


(cont.) 


ROCKET SAMPLE OUTPUT 


F WEIGUT WDOTP 


358. 
362. 


362. 


362. 
362. 
362. 
362. 


362. 


362. 
362. 
362. 
362. 


362. 
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51.94 


INS WT 


29. 


30 


. WAS 


83 


met NW @ & UT UT OF ~~) oO Ww 


153. 


1110.933 
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Table 3: (cont.) 
ROCKET SAMPLE OUTPUT 


TIME eC PE F WEIGHT CF WDOT EX ISP 
004 155.57 3.91 1603.81 -12 1.191 02 155.91 
00s 167.45 3.85 1746.08 13° 1.204 -02 358.41 
00s 191.48 3.75 2036.21 15 1.228 -03 162.96 
006 215.74 3.66 2353.98 18 1.260 .04 167.19 
007 240.11 3.59 2676.77 20 1.288 OS 371.24 
007 264.49 3.52 3002.53 23 3.311 -06 175.05 
008 288.80 3.46 3329.74 25 1.332 .07 178.621 
009 312.95 3.56 3674.16 28 1.356 -08 181.99 
009 336.89 3.83 4033.64 31 1.383 .09 1865.31 
010 360.55 4.10 4389.03 33° 1.406 -10 188.54 
010 383.89 4.37 4739.64 36 1.426 -12 191.66 
O11 406.88 4.63 5084.88 39 1.443 -13 194.63 
012 451.65 5.14 5757.28 45 1.472 -16 200.10 
014 494.66 5.63 6403.28 51 1.495 -20 204.95 
015 535.78 6.10 7020.91 57 1.514 -23 209.24 
016 $74.94 6.54 7609.00 63 1.529 -27 213.03 
018 612.09 6.97 8167.05 70 «1.541 -31 216.40 
019 647.24 7.37 8695.02 76 1.552 .36 219.40 
020 680.42 7.74 9193.27 83 1.561 -41 222.08 
021 711.65 8.10 9662.43 90 1.568 -46 224.49 
023 JAL.O1 8.43 10103.34 -97 1.575 -52 226.66 
024 768.55 8.75 10517.02 1.04 1.581 -56 228.62 
025 794.35 9.04 10904.56 1.11 1.586 -62 230.40 
026 818.49 9.31 11267.14 1.18 1.590 .67 232.02 
029 862.10 9.81 11922.16 1.32 1.597 79 234.86 
032 900.05 10.24 12492.06 1.47 1.603 91 237.25 
034 932.96 10.62 12986.34 1.62 1.608 1.04 239.29 
037 961.43 10.94 13414.09 1.76 1.611 1.18 241.05 
039 986.03 11.22 13783.48 1.92 1.615 1.31 242.57 
042 1007.24 11.46 14102.08 2.07 1.617 1.45 243.90 
044 1025.51 11.67 14376.50 2.22 1.619 1.59 245.07 
047 «1041.23 11.85 14612.75 2.37 1.621 1.74 246.10 
0S0 1054.76 12.00 14815.98 2.53 1.622 1.88 247.02 
052 1066.40 12.14 14990.82 2.68 1.624 2.03 247.84 
OSS 1076.41 12.25 15141.17 2.84 1.625 2.18 248.57 
060 1092.43 12.43 15381.987 3.15 1.626 2.48 249.83 
065 1104.32 12.57 15560.39 3.46 1.627 2.79 250.87 
073 1116.71 12.71 15746.61 3.94 1.629 3.25 252.12 
083 1126.81 12.82 15898.38 4.57 1.630 3.87 253.37 
106 1137.17 12.94 16054.16 €.090 4 633 §.29 255.19 
183 1148.30 13.07 16222.20 10.79 1.632 10.05 257.68 
249 1155.79 13.16 16335.24 14.97 1.632 14.21 258.55 
255 1140.96 12.99 16112.51 18.33 1.631 14.57 258.60 
263° 1128.01 12.84 15918.07 1$.80 1.630 15.05 258 .66 
273° #1117.77 12.72 15764.40 16.42 1.629 15.67 258.71 
311 1107.56 12.61 15611.44 18.74 1.628 17.99 258.83 
322 1107.41 12.61 15609.22 19.35 1.628 18.60 258.86 
532 1118.03 12.73 18770.64 32.10 1.629 33.29 259.19 
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TIME 


ANDAWUUTUU UD ee W WD pe pe 


NAAHANANARAAADADADANADADAAAD 


736 
-982 
-228 
576 
.006 
-539 
-€00 
-112 
645 
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Table 3: 


F 
18927. 
16081. 
16235. 
16393. 
16553. 
16649. 
16652. 
16627. 
16498. 
16334. 
16169. 
15996. 
15818. 
15669. 
15503. 
1$351. 
15200. 
15045. 
14865. 
11683. 
14502. 
34335. 
14149. 
13904. 
13664. 
13427. 
12964. 
12517. 
12083. 
11258. 
10484. 

9080. 
7848. 
6765. 
5815. 
4980. 
A248. 
3604. 
3085. 
2644. 
2261. 
1936. 
1667. 
1433. 
1230. 
1054. 

902. 

769. 
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(cont.) 
ROCKET SAMPLE OUTPUT 


WEIGHT 
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WDOT EX 
-77 


260. 


260. 


Table 3: (cont.) 
ROCKET SAMPLE OUTPUT 


TIME PC PE E: WEIGHT CF WDOT EX ISP 

6.463 73.46 54 655.25 393.14 1.026 392.88 260.02 

6.473 56.63 79 470.06 393.14 -955 392.92 260.01 

6.484 43.66 : . -874 392.95 260.00 

6.494 33.65 31 227.84 393.14 .779 392.97 259 .99 

6.509 22.78 74 119.40 393.14 -603 392.99 259.99 
AVERAGE THRUST FROM -000 SEC. TO 6.509 SEC. WAS 15696.570 Ibf 
TOTAL IMPULSE= 102172.50 SEC. 

THRUST VS. TIME 
THRUST 
(1000 Ibf) 
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Table 3: (cont.) 
ROCKET SAMPLE OUTPUT 


FRESSURF VS. TIME 


PRESSURE 
(1000 psi) 


-_ ete ©» ee 
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IV. RESULTS AND CONCLUSIONS 


A. RESULTS 

In order to verify the accuracy of the modifications 
made to NEWPEP, several problems were run using the old 
version, the updated version and the NASA SP-273 code [Ref. 
6]. The results for the following two problems are 
summarized in Table 4. 

-1. A mixture of 6% by weight H, and 94% by weight O, 
was burned at a chamber pressure of 3.447 MPa (500 psi) and 
expanded to an ambient pressure of 0.1014 MPa (14.7 psi). 

2. A mixture of 85% by weight ammonium perchlorate and 
15% HTPB (Sinclair) was burned at 6.8947 MPa (1000 psi) and 
expanded to an ambient pressure of 0.1014 MPa (14.7 psi). 
There were no significant differences in any of the output 


values. 


B. CONCLUSIONS 
PEP93 is an interactive program that works well ona 


personal computer. It allows the user to evaluate a wider 


variety of performance/design problems while using an 


existing, familiar program. 
ROCKET was successfully modified to allow the use of 
multiple grains. This greatly enhances the ability to 


predict the performance of actual rocket motors. 





Table 4: 
SAMPLE COMPARISON OF NEWPEP/PEP93/NASA SP-273 


__|AP 85% by weight ___ 
{_ [HTP 15% by weight_ 
_,NEWPEP |PEPS3___|NASA 


2344) 2044) 2644.19 
8} 2192.79) 2199.18 
10.45} 10.146 





___1s01] 1501 
8.091} 8.0894 
__ 1575] 1.575 
___ 24a) 243.8 
1513.5) ‘1514 
___8.361] | 8.3596 


~ 9.105] 8.104 


157M] 1.571 
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APPENDIZ A 


USER'S GUIDE FOR PEP93 


(Changes in Bold) 


A. GENERAL 

1. The files NEWPEP.FOR, NEWEQUIL.FOR and NEWDSIGN. FOR 
are needed to create NEWPEP.EXE. In addition, several 
common blocks are needed. BLANK.INC, CHARA.INC, IBRIUM.INC, 
SCRATCH. INC and SIZE.INC. 


2. The data files JANNAF.DAT and PEPCODED.DAT are 


needed to run PEP93. 


B. INPUT 
a bee Up to 15 ingredients may be input by the user. 
These ingredients may be coded ingredients from the 
PEPCODED.DAT file or user defined ingredients, or both. 
2. The user is prompted for the following inputs when 
PEP93 is run: 
File name (up to 30 characters) 
b. Run title (up to 17 characters) 
ec. Vitiated air option 
(1) Number of coded ingredients for vitiator 
(2) Number of user defined ingredients for 
vitiator 


(3) Ingredient code numbers from PEPCODED.DAT 
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User defined ingredient name, composition, 
heat of formation and density 

(5) Chamber pressure 

(6) Chamber inlet temperature 

(7) Ingredient weights 


Number of coded ingredients for combustor 


Number of user defined ingredients for combustor 


Number of runs 

Density exponent used in calculating density 

(defaults to 1.0) 

Propellant temperature to adjust the system 

enthalpy for heating or cooling 

(default is 298K) 

Quadratic coefficients for solid specific heat. 

Used to adjust the system enthalpy for heating 

or cooling (defaults to 0.3 cal/gm-K) 

Nozzle discharge coefficient (defaults to 1.0) 

Options 

(1) Delete exit calculations 

(2) Include ionic species in calculations 

(3) Include boost velocities and nozzle design 
data 
Input pressures in psi instead of MPa 
Increase numerical precision of species 


concentrations 





(6) Output a list of all combustion species 
considered 
(7) Fix temperature 
(8) Debug options 
(a) Thermo data at every guess 
(b) Values of J, M, VF, PR, VA in 
subroutine TWITCH 
(c) Species composition every guess 
(d) Log of equilibrium constant every 
guess 
(e) Classification of species each 
iteration from TWITCH 
(9) Exact or approximate throat calculations 
(10) Normalize or not normalize to 100 grams 
(11) Calculate non-ideal expansion performance 


(12) Calculate chamber Mach number 


1. Ingredient code numbers from PEPCODED.DAT 

m. User defined ingredient name, composition, heat 
of formation and density 

n. Chamber pressure (MPa) 

o. Exhaust pressure or area ratio and ambient 
pressure (see note 2) 

p. Ingredient weights 

q. Ratio of chamber area to throat area 

Notes: 

a. Format for user defined ingredients 
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(1) Ingredient name - A30 


(2) #atoms of element/symbol of element - 6(I3, 
A2) 
(3) Heat of formation(cal/gram) - F5.0 
(4) Density(gm/em?) - F6.0 
b. If the non-ideal expansion option is not chosen, 
the exhaust pressure is used for input "o" and is also equal 
to the ambient pressure. If the non-ideal expansion option 


is chosen, the ratio (exit area/throat area) and the ambient 


pressure are input at "o". 
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C. SAMPLE INPUT 


TYPE IN OUTPUT FILE NAME (DEFAULT: PEPOUT.DAT) 

If you wish to direct to a printer, type in 
printer name such as "LPT1" for the IBM PC 

or "LPAO:" for the VAX. 

To send output to screen, type "TT;" (VAX & IBM PC) 


TYPE RUN TITLE (char*17) OR type quit. 

sample 

DO YOU WANT TO CALCULATE COMPOSITION AND ENTHALPY OF 
VITIATED AIR? 1=YES,2=NO 

1 

INPUT NO. OF CODED ING., NO. OF USER DEFINED ING. FOR 
VITIATOR 

2,0 


READ IN 2 INGREDIENT CODE NUMBERS SEP. BY COMMAS 
44,729 


READ IN CHAMBER PRESSURE (MPa), CHAMBER INLET TEMPERATURE 
(K), wt1,wt2, +etc. 2.068,700.,96.,4. 


BEGIN ENGINE DATA INPUT. 


NO. OF CODED ING., NO. OF USER DEFINED ING., AND NO.OF RUNS 


NOTE: INPUT VITIATED AIR AS A USER DEFINED INGREDIENT. 
1,1,1 : 
DENSITY EXPONENT (DEFAULTS TO 1.0) = 


PROPELLANT TEMPERATURE (DEFAULTS TO 298.0) = 


ENTER THREE VALUES FOR CSUBP FIT 
FIRST VALUE = 0.0 DEFAULTS TO 0.30 


NOZZLE DISCHARGE COEFFICIENT (DEFAULTS TO 1.0) = 


OPTIONS; 1-DELETE EXIT CALCULATIONS 

2-INCLUDE IONIC SPECIES IN CALCULATIONS 

3-INCLUDE BOOST VELOCITIES AND NOZZLE DESIGN DATA 

4-INPUT PRESSURES IN PSI INSTEAD OF MPA 

5-"N” ORDERS OF MAGNITUDE SPECIES PRECISION 
INCREASE 

6-OUTPUT A LIST OF ALL COMBUSTION SPECIES 
CONSIDERED 

7-DEBUG OPTIONS 

8-"P-T-H-S" MAP OPTION 

9-"O0"=sexact throat "1"=approx throat "2"=exact 
with throat composit. 

0-"1"=do not normalize weights to 100 gms 
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"0"=normalize 
A-CALCULATE NON-IDEAL EXPANSION PERFORMANCE 
B-CALCULATE CHAMPER MACH NUMBER 
1234567890AB**put O-no or 1-yes under number (opt 5 takes 0 
or"™n") 
000000000011 


READ IN 1 INGREDIENT CODE NUMBERS SEP. BY COMMAS 
714 


VITIATED AIR IS AUTOMATICALLY ENTERED AS FIRST USER DEFINED 
INGREDIENT 


READ IN CHAMBER PRESSURE (MPa) ,AREA RATIO, AMBIENT PRESSURE 
(MPa), WT1,WT2, +ETC. 

INCLUDE DECIMAL POINT AND SEPARATE BY COMMAS 

INPUT A ZERO CHAMBER PRESSURE TO ADD INGREDIENTS AGAIN 
INPUT A -1.0 TO END PROGRAM 

1.379,5.,14.7,7.,93. 


INPUT RATIO OF CHAMBER AREA TO THROAT AREA TO. 
CALCULATE MACH NUMBER AT STATION 4 
1.5 


TYPE RUN TITLE (char*17) OR type quit 


quit 
Stop - Program terminated 


D. SAMPLE OUTPUT 


See Table 2. 





APPENDIX B 


USER’S GUIDE FOR ROCKET 


A. INPUT 
1. The data file CARD.DAT is needed to run ROCKET and 
consists of the following variables. 
Title 
Ath 
Pi 


ETH 


pressure 1 erosion rate at pressure 1 


pressure 2 erosion rate at pressure 2 
etc. etc. 
G# 


° 
Yr 


where 
Ath Throat area (in?) 
Enter a positive area if it is known 
Enter a negative area for a first guess if unknown 


Enter zero if you have no first guess (very slow) 


Ae Nozzle exit area (in?) 





Vol 


Not 


Pd 


td 


ET# 


Total motor volume without propellant or igniter (in?) 
Ambient propellant temperature (°F) 
Ambient pressure (psia) 
Throat plug closure blowout pressure (psia) {typically 
35} 
Initial motor pressure (psia) 
Nozzle ioss coefficient 
Generally use y, for shifting equilibrium flow 
(PEP93) 
Design pressure for sizing the throat (psia) 
Throat erosion delay time (sec) 
Positive value: fixed rate from t=0 (i.e. 0.0005 
in/sec) 
Negative : equals the time delay for onset of erosion 
Negative or 0: must include pressure vs. erosion rate 
table (Erosion rates are radial in in./sec) 
Number of points for throat erosion table 
Number of grains 
Propellant burning rate at P,,, (in./sec) 
Burning rate exponent (r=aP,*) 
Propellant temperature sensitivity (%/°F) 
Reference pressure for the specified burning rate 
(psia) 
Characteristic exhaust velocity at P,., (ft/sec) 


Pressure correction exponent for C* 
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Density of grain(lbm/in?) {propellant,igniter or 
insulation} 

Ignition time (sec) {0 if no igniter} 
Ignition time delay (sec) {0 if no igniter} 
Number of points in web burned vs. burning surface 
area table 

Web cistance burned (in) 

Propellant burning surface area (in?) 

Exposed insulation area (in?) 

Identification number for grain 

1:igniter 

2:insulation 

3:propellant 

Insulation ablation rate 


2. Notes: 

a. All input data is formatted F12.6 except: 
(1) Title - no format required 
(2) ET#,GH# - format I3 

b. Enter 0 in columns 3,6 and 9 of the last row. 

c. Igniters: 
(1) © always equal to 1 
(2) pts is # of points in time vs. m curve 
(3) Wb=time (sec) 
(4) Ab= m (lbm/sec) 

d. Insulation: 


(1) r,n,I, and P.,.- refer to the propellant 
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grain, not to the insulation. They are used to calculate a 
reference web which can be used to describe the amount of 
insulation exposed during a firing. C*,X,p, ,t, and At, are 
for the insulation. 

(2) Wb is the web of the associated propellant 


grain. Ab is the exposed insulation area. 


B. OUTPUT 
1. Input data is printed with all values labeled. 
2. Output data consists of the following: 

a. Total propellant weight and igniter weight. 

b. A table of time, chamber pressure, change in 
pressure, total motor volume, propellant weight consumed, 
rate of mass addition due to propellant burning, insulation 
weight consumed, rate of mass flow through the nozzle and 
throat area. 

c. Average and maximum pressure for run. 

d. A table of time, chamber pressure, exit 
pressure, thrust, coefficient of thrust, rate of mass flow 
through the nozzle and specific impulse. 

e. Average thrust and total impulse. 

£. Plots of thrust versus time and pressure versus 

time. 
3. Time, pressure, and thrust are also put into the 


data file PLOT.DAT 
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C. SAMPLE INPUT 


See Table 5. 


D. SAMPLE OUTPUT 


See Table 3. 
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15percent 


3 


11 


11 


0 


8.657700 
14.700000 


0.580000 
0.062300 
3 

-000000 
- 380900 
- 761800 
- 142800 
-523700 
-904600 
- 285500 
- 6664100 
-047400 
- 428300 
- 809200 
- 580000 
-062300 


-000000 
- 380900 
- 761800 
-142800 
-523700 
-904600 
- 285500 
-666400 
-047400 
- 428300 
- 809200 
- 000000 
-962300 
1 
0.000000 
9.250000 
0 0 


DM WWWNNE HH OCS FPOWVUWUWNNH Se OOO 


Table 5: 
ROCKET SAMPLE INPUT 


105.450000 8042.477000 
0.960600 


0.350000 
9.000000 


1448. 
1519. 
1542. 
1557. 
1563 
1562. 
1552. 
1534. 
1506. 
1473. 
1431 


300000 
500000 
500000 
300000 


-900000 


300000 
500000 
500000 
200000 
800000 


- 200000 
- 350000 
-000000 


. 300000 
-500000 
-500000 
- 300000 
-900000 
- 300000 
-500000 
-500000 
- 200000 
- 800000 
- 200000 
-000000 
-000000 


-400000 
- 900000 


1.144300 


0.200000 
0.000000 


0.200000 
0.000000 


0.0nn00n 


70.000000 10.130000 
1030.900000 0.000500 


1030.9000000 5106.660000 


1030.000000 S106.660000 


0.050090 


0.000000 5106.660000 


47 


35.000000 


0.045000 


0.045000 


9.045000 


APPENDIX C 


PEP93 SUBROUTINE EQUIVALENCE 


SUBROUTINE EQUIVALENCE (equiv) 
INCLUDE ‘'size.inc’ 
INCLUDE ‘blank.inc’ 
INCLUDE ‘chara.inc’ 
INCLUDE ’scratch.inc’ 
INCLUDE ‘ibrium.inc’ 
real ostoic,oin, equiv 
common/llist/atwt (100) 
ostoic=0.0 
do 10 j=i1,is 


if (aspec(j) .eq. ‘'Ag’) ostoic=ostoic+alp(j) 
if(aspec(j) .eq. ’Al’) ostoic=ostoic+alp(j)*3.0/2 
if (aspec(j) .eq. 'B’ ) ostoic=ostoic+alp(j)*3.0/2 
if (aspec(j) .eq. ‘'Ba’) ostoic=ostoic+alp(j) *2.0 
if (aspec(j) .eq. 'Be’) ostoic=ostoic+alp(j) 
if(aspec(j) -eq. 'Bi’) ostoic=ostoic+alp(j)*3.0/2 
if(aspec(j) .eq. ‘'C’ ) ostoic=ostoic+alp(j) *2.0 
if (aspec(j) .eq. ‘'Ce’) ostoic=ostoic+alp(j) *2.0 
if(aspec(j) .eq. ‘’Cr’) ostoic=ostoic+alp(j) *3.0 
if (aspec(j) .eq. 'Cs’) ostoic=ostoic+alp(j)*3.0/2. 
if(aspec(j) .eq. ‘Cu’) ostoic=ostoic+alp(j) 

if (aspec(j) .eq. ‘'Fe’) ostoic=ostoic+alp(j) *3.0/2 
if (aspec(j) .eq. ‘H’ ) ostoic=ostoic+alp(j)/2.0 
if(aspec(j) .eq. 'Hg’) ostoic=ostoic+alp(j) 
if(aspec(j) .eq. ‘'K’ ) ostoic=ostoic+alp(j) *2.0 
if(aspec(j) .eq. ‘Li’) ostoic=ostoic+alp(j) 

if (aspec(j) .eq. 'Mg') ostoic=ostoic+alp(j)*2.0 
if(aspec(j) .eq. ‘Mn’) ostoic=ostoic+alp(j)*2.0 
if (aspec(j) .eq. 'Mo’) ostoic=ostoic+alp(j)*3.0 
if(aspec(j) .eq. ‘'Na’) ostoic=ostoic+alp(j)*2.0 
if (aspec(j) .eq. ‘Ni’) ostoic=ostoic+alp(j) 

if (aspec(j) .eq. ‘'P’' ) ostoic=ostoic+alp(j)*5.0/2. 
if(aspec(j) .eq. ‘Pb’) ostoic=ostoic+alp(j)*2.0 
if(aspec(j) .eq. ‘Rb’) ostoic=ostoic+alp(j)*2.0 
if(aspec(j) .eq. 'S' ) ostoic=ostoic+alp(j)*2.0 
if(aspec(j) .eq. ‘'Si’') ostoic=ostoic+alp(j)*2.0 
if (aspec(j) .eq. ‘'Sn’) ostoic=ostoic+ralp(j) *2.0 
if(aspec(j) .eq. ‘'Th’) ostoic=ostoic+alp(j)*2.0 
if (aspec(j) .eq. ‘’Ti') ostoic=ostoic+alp(j)*2.0 
if (aspec(j) .eq. ‘U’ ) ostoic=ostoic+alp(j)*2.0 
if(aspec(j) .eq. ’V’ ) ostoic=ostoic+alp(j)*5.0/2. 
if(aspec(j) .eq. ‘'W’ ) ostoic=ostoic+alp(j)*3.0 
if (aspec(j) .eq. ‘'Zn’) ostoiceostoic+alp(j) 
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if (aspec(j) .eq. ‘Zr’) ostoic=ostoic+alp(j)*2.0 
if(aspec(j) .eq. '0’)oin=alp(j) 
10 continue 
equiv=ostoic/oin 
return 
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